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ABSTRACT 


The  Istituto  Elettrotecnico  Nazionale  time  scale  UTC(IEN)  is  generated  by  a 
commercial  cesium  clock  selected  among  the  five  presently  at  disposal  and  is 
maintained  in  agreement  within  5  microseconds  with  UTC  by  means  of  a  phase  mi¬ 
crostepper  . 

The  problem  of  the  influence  of  the  environmental  parameters  on  atomic  clocks  has 
been  investigated  and  recent  improvements  in  the  international  synchronization 
systems  have  shown  more  clearly  the  presence  of  seasonal  effects  on  clocks. 

On  the  IEN  cesium  clock  ensemble,  long  term  frequency  instabilities  and  seasonal 
variations  have  been  detected  and  a  relationship  between  humidity  and  frequency 
variations  has  been  analyzed.  This  investigation  has  been  also  extended  to  the 
cesium  clocks  kept  in  the  (ISPT)  Laboratory  in  Roma. 

An  evaluation  of  the  annual  repetition  of  these  effects  has  been  performed  by 
means  of  statistical  processes  to  settle  a  better  use  of  microstepper  corrections 
and  optimize  the  long  term  stability  of  UTC(IEN)  time  scale. 


1 .  Introduction 

The  influence  of  environmental  conditions  on  the  commercial  cesium  clock 
behaviour  has  been  analyzed  in  the  past  and  reported  in  the  literature  [ 1,2, 3, 4] 
where  sensitivity  to  temperature  and  humidity  variations  was  pointed  out. 

In  spite  of  the  careful  investigations  already  performed,  it  has  not  yet  been 
possible  to  relate  clearly  the  amount  of  the  rate  changes  of  the  clocks 
considered  to  the  variations  introduced  in  their  environmental  conditions.  More¬ 
over,  the  results  we  have  at  disposal  [2,5],  have  demonstrated  that  each  clock 
must  be  analyzed  separately  to  find  out  a  "best  estimate"  of  its  long  term  beha¬ 
viour  where  sensible  changes  in  its  environment  can  be  expected.  This  is  in¬ 
strumental  in  maintaining  a  time  scale  of  a  laboratory  in  close  agreement  with 
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the  international  atomic  time  scale  if  the  source  of  the  local  UTC  is  a  cesium 
clock,  whose  rate  is  periodically  adjusted,  and/or  if  it  is  a  paper  time  scale 
computed  from  a  clock  ensemble. 


2 .  Long  term  behaviour  of  the  UTC (TEN)  time  scale 

The  Istituto  Elettrotecnico  Nazionale  time  scale  UTC(IEN)  is  realized  by  means  of 
five  commercial  cesium  beam  standards  [6],  that  are  kept  at  a  temperature  of 
24+0. 5° C  in  a  room  located  12  meters  underground. 

The  laboratory  master  clock  is  chosen  according  to  its  long  term  stability  beha¬ 
viour  as  deduced  by  the  IEN  intercomparisons  and  from  the  bi-monthly  rates  versus 
the  international  atomic  time  scale  data  published  by  the  BIPM[7].The  master  clock 
frequency  is  supplied  to  a  phase  microstepper  whose  rate  is  periodically  adjusted 
to  maintain  UTC (IEN)  in  agreement  within  5  us  with  UTC  established  by  the  BIPM. 

Since  mid-1985,  the  IEN  time  scale  is  compared  with  UTC  by  means  of  the  Global 
Positioning  System  satellites,  according  to  the  common  view  measurement  schedule. 
This  change  in  the  international  synchronization  link  used  by  IEN  as  well  as  by 
many  other  laboratories  has  significantly  improved  the  accuracy  of  the  time 
comparisons  [8]  allowing  new  investigations  on  periodic  rate  changes  in  atomic 
clocks . 

In  Fig.  1,  the  equipment  set  up  for  the  IEN  time  scale  realization  is  depitected, 
meanwhile,  in  Fig.  2,  its  behaviour  against  UTC  time  scale  in  the  last  three 
years  is  shown. 

Table  1  reports  the  IEN  master  clock  rate  corrections  performed  by  means  of  the 
phase  microstepper  to  maintain  UTC(IEN)  in  close  agreement  with  UTC  in  the  period 
1985-1987. 

Table  1 


MJD 

Date 

Master  Clock 

UTC (IEN)  -  Master  Clock 
(ns/d) 

46213 

1985 

May 

28 

Cs5061A-004  sn893 

+  34 

46473 

1986 

Feb 

12 

II 

+  10 

46723 

1986 

Oct 

20 

II 

+32 

46913 

1987 

Apr 

28 

II 

+  16 

46983 

1987 

Jul 

8 

II 

+  28 

Looking  carefully  at  the  shape  of  the  UTC(IEN)  graph  it  can  be  seen  that  a 
seasonal  effect  is  present  and  that  the  rate  compensation  performed  on  the  phase 
microstepper  has  not  eliminated  this  effect. 


Moreover  the  time  scale  generated  directly  by  the  master  clock,  as  reported  in 
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Fig.  3,  would  have  been  more  uniform  without  the  use  of  the  microstepper,  but 
with  a  total  time  error  exceeding  our  stated  limits.  This  consideration  showed 
the  importance  of  a  careful  characterization  of  the  long  term  behaviour  of  every 
IEN  clock  to  find  out  if  seasonal  effects  were  detectable  in  order  to  optimize 
the  time  scale  corrections. 

The  environmental  conditions  of  the  clock  room,  as  regards  temperature  and  ab¬ 
solute  humidity,  measured  in  the  period  1985-1987  are  reported  in  Fig.  4.  It  can 
be  noticed  that,  while  the  temperature  is  stable  within  +0.5°C,  the  absolute  hum¬ 
idity  can  vary  from  about  8  gm  to  18  gm  (about  40%  to  80%  relative  humidity) 
(*). 

This  is  due  to  the  fact  that  in  this  room  only  the  temperature  is  controlled  by 
means  of  a  heater  switched  on  automatically  by  a  temperature  sensor  meanwhile  no 
humidity  control  is  made  and  therefore  the  external  variation  is  followed  with 
nearly  a  10  days  delay. 

A  first  attempt  in  analyzing  the  long  term  behaviour  against  UTC  of  some  cesium 
clocks  at  IEN  in  Torino  and  at  the  Istituto  Superiore  Poste  e  Telecomunicazioni 
(ISPT)  in  Roma  [9]  showed  the  presence  of  seasonal  effects.  However  it  turned  out 
necessary  to  improve  the  investigation  tools  by  adopting  suitable  mathematical 
models  and  statistical  treatments. 


3 .  Measurements  set-up  and  synchronization  links  used 

The  time  differences  between  UTC(IEN)  and  the  IEN  cesium  clocks  analyzed  and 
reported  in  this  paper  are  ten  day  averages  of  the  values  measured  daily  at  OOh 
UTC  with  a  1  ns  resolution. 

Concerning  the  ISPT  cesium  clocks,  the  measurement  resolution  is  the  same,  but 
they  are  referred  to  UTC(IEN)  by  means  of  daily  television  measurements  having  a 
typical  standard  deviation  of  50  ns  (l<r).  This  laboratory  has  both  humidity  and 
temperature  controls. 

As  the  IEN  has  only  commercial  cesium  standards,  all  the  clocks  readings  have  be¬ 
en  referred  to  UTC,  by  means  of  the  corrections  of  UTC(IEN)  published  by  BIPM, 
assuming  that  the  international  atomic  time  scale  is  less  sensitive  to  seasonal 
effects  being  made  up  by  a  large  number  of  clocks  kept  in  different  locations  and 
environments . 

The  IEN  and  ISPT  temperature  and  absolute  humidity  data,  have  also  been  averaged 


(*)  In  practice,  in  the  clocks  room  the  relative  humidity  in  %  is  measured.  These 
values  have  been  converted  into  absolute  humidity  data,  expressed  in  gm  , 
using  the  conversion  coefficients  given  by  the  Smithsonian  Physical  tables, 
related  to  the  room  temperature. 
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over  ten  days  from  the  daily  readings  taken  by  the  data  acquisition  systems  with 
resolution  of  . 1  K  and  of  5%. 


4 .  Mathematical  models  for  seasonal  variations 

When  dealing  with  seasonal  variations  it  is  helpful  to  use  a  simple  but  adequate 
mathematical  model  to  extract  quasi-periodic  variations  from  drifts  and  statisti¬ 
cal  fluctuations.  In  this  paper  two  different  attempts  have  been  performed. 

With  the  first  approach  we  approximate  the  long  term  behaviour  of  the  time 
difference  between  a  time  reference  and  a  cesium  clock  with  the  function 


-  A  4-  bpc  +  —  C  -  d  CO-S  uux.  (1) 

where:  x  is  expressed  in  days,  f(x)  in  microseconds ,tu  ,  the  angular  frequency,  is 
equal  to  2ff/(365  days).  The  variable  x  ranges  from  0  to  k*365  while  uJx  varies 
from  0  to  k*2tr,  where  k  is  an  integer  number.  The  parameters  a,b,c,d  can  be 
calculated  by  subsequent  numerical  derivations,  assuming  that  the  period  of  the 
sinusoidal  part  is  exactly  365  days  and  analyzing  the  data  relative  to  an  integer 
number  of  periods  to  be  sure  that  the  mean  value  of  cosu»x  is  null. 

This  method  gives  satisfactory  results  but  the  differentiation  magnifies  the 
irregularities  of  the  data  and  therefore  it  is  necessary  a  smoothing  process  at 
each  derivation  step.  On  the  contrary,  a  method  based  on  an  integration  process 
does  not  present  this  kind  of  problem. 

Considering  always  an  integer  number  of  periods,  we  can  calculate  the  first  terms 
of  a  Fourier  series,  obtaining: 


cliuyc  =  2*3 
S’nuoe  d  ujuc  = 
^-Gc)  cnsuoc  dulx  = 

J  T?Gt)  C03  du JX  • 


2.tr 


2-T 


3  >juz 

— 

uU 

—  ol  IT  ■+ 


TTC 


2,  TTC. 


uu 


2.  uu‘ 


(2) 


Operating  by  numerical  integration  and  substitution  in  Equations  (2),  the  values 
of  the  coefficients  a,b,c,d  are  obtained. 

It  turns  out  that  these  two  methods  give  almost  the  same  results. 

The  experimental  data  can  also  be  fitted  to  an  arbitrary  function,  for  example 
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with  the  Gauss-Newton  method  [ 10]  which  satisfies  the  least-squares  approximati¬ 
on.  This  iterative  process  has  been  applied  to  the  experimental  data  and  stopped 
when  the  relative  parameter  variations  were  lower  than  10  .  It  has  been 
estimated  more  suitable  to  use  in  the  calculations  the  normalized  frequency 
departures  data  of  each  clock  considered.  Therefore  the  fitting  curve  defined  in 
(1)  is  modified  into 


where  the  numerical  factor  allows  to  express  h(x)  in  10  .  The  parameters  values 
so  obtained  are  the  best  approximation  as  far  as  the  least-squares  regression  is 
concerned.  Moreover  the  possible  initial  phase  shift  has  been  introduced  and 
evaluated  to  avoid  any  subjective  estimate  of  the  time  origin  of  the  sine 
function.  Consequently  the  final  fitting  function  is 


J 


uu  _5<n  uu 


(dc  +  q)  +  h> 


where  e  is  expressed  in  days. 


(4) 


As  a  last  step,  the  statistical  fluctuations  can  be  computed  as  the  residual  bet¬ 
ween  the  experimental  data  and  the  approximating  curve  h(x) .  These  fluctuations 
can  be  analyzed  by  means  of  the  Allan  variance  to  characterize  the  kind  of  fre¬ 
quency  noises  affecting  the  clocks.  An  example  of  this  statistical  computation, 
performed  on  Cs9,  is  given  in  Table  2.  A  white  frequency  noise  can  be  evinced 
from  this  data. 


Table  2 


cry  (r)  [10  14  ] 

10 

5.8 

20 

6.1 

40 

4.8 

80 

3.1 

160 

2.4 

The  amount  of  the  seasonal  fluctuations  has  been  determined  as  well  through  a 
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second  approach  based  only  on  statistical  treatments.  A  cesium  standard  phase 
variation  with  respect  to  a  reference  is  assumed  to  be  of  the  form 


n*._, 

ic! 


+£  a.  It)  si n  L ujmt  +  6jt) + Qj 


+ 


(5) 


where  oJ0  and  ujm  are  angular  frequencies  of  operation  and  of  periodic  seasonal 
effects  respectively ,  -O^k !  are  long  term  phase  drift  coefficients ,  Q^t)  ,  Q^Ct) 
and^p(t)  are  random  variables.  Seasonal  variations  can  be  evaluated  by  computing 
variances  if  the  polynomial  drift  is  unimportant  either  because  it  is  very  small, 
in  the  sample  time  of  interest,  compared  with  periodic  modulations,  or  it  is 
removed  by  using  Hadamard  variances  on  sets  of  a  suitable  number  M  of  adjacent 
samples  weighted  with  binomial  coefficients  [11,12]. 

In  the  following,  the  expressions  concerning  the  simple  special  cases  of  the 

two-sample  variance  O'1  and  of  the  Hadamard  variance  <crt  (M  t  with  M  =  3 

7  "BC  7  y 

are  reported. 

For  sake  of  simplicity  X  cc  (t)  is  substituted  with  a(t)  whose  mean  value  and 
variance  are  and  <r£  respectively,  moreover  the  correlation  coefficients  of 
the  fluctuations  of  a(t)  and  0(t)  are  supposed  to  be  zero  for  the  time  intervals 
of  the  order  and/or  longer  than  half  a  year.  The  same  hypothesis  is  assumed  for 
all  the  other  covariance  coefficients. 

For  the  Allan  variance  case,  if  contributions  from  all  n  ^are  negligible,  the 
following  expressions  are  obtained: 


■h  C052  90  [  \( 


if  vjJ [r  =  Tf  ,  and 


-h  & 


-O', 


& 


(  2.  cclsu^-c  +  '/^  eczs  2,  uJMr)|(6b) 


if  TT  • 

In  (6a)  and  (6b)  O'g  is  the  variance  of  0(t)  and  <3^  is  the  Allan  variance  of 
the  process  (iXt). 
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For  the  Hadamard  variance  (M  =  3),  under  the  same  assumptions  stated  before,  the 
following  expressions  hold: 


(7a) 


(7b) 


if  o>  TT  ■ 

Useful  simpler  expressions  can  be  obtained  to  evaluate  the  effects  of  small 
fluctuations  in  a(t)  and  9(t),  that  is  when  %  «  1  and  oi  //<*  cc  1 .  If  these 
fluctuations  are  entirely  negligible,  the  substitutions  C*  =  c?g  =  0  in  (6)  and 
(7)  lead  to  the  following  results  for  periodic  contributions: 


if  =  Ti 


if  4-  T 


i  f  £  =  TT 


if  ±  r r 


(8  a) 

(8b) 

(9a) 

(9b) 


In  the  special  case  iUm"C  -  ''T  the  two  variances  exhibit  the  same  dependence  on  Q0  , 
whereas  for  £yt  these  simples  expressions  show  that  the  Hadamard  variance  is 
more  selective.  The  curves  of  and  of  (0(1  Jit r)  >  / [32  ( A,  / LUcc  )*] 
referring  to  a  sinusoidal  deterministic  process  are  compared  in  Fig.  5. 
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5.  Evaluations  of  the  seasonal  variations 


The  mathematical  and  statistical  processes  mentioned  in  Section  4  have  been 
applied  to  the  cesium  clocks  kept  at  IEN  and  at  ISPT.  In  Table  3  the  model,  the 
serial  number  and  the  period  of  observation,  in  Modified  Julian  Date  (MJD),  are 
reported  for  each  clock. 

Table  3 


Laboratory  Clock  Manufacturer  Serial  number  Abbreviation  Observation 


and  model 

used 

period 

from-to 

(days)  (MJD) 

ISPT 

HP5061A/004 

597 

Csl 

450 

46229-46679 

HP5061A/004 

1815 

Cs3 

570 

46489-47059 

IEN 

HP5061A 

609 

Cs5 

660 

46229-46889 

HP5061A/004 

893 

Cs6 

730 

46229-46959 

0SA3200 

84 

Cs9 

830 

46229-47059 

The 

normalized  frequency 

departure  versus  UTC  with 

a  constant  value  subtracted 

and 

the  frequency  drift  removed  are  reported  for  each  clock  in  Fig. 

6  and  Fig.  7. 

The 

results  obtained  by 

fitting  the 

experimental 

data  using  the 

Gauss-Newton 

method  have  been  plotted  over  the  same  graphs,  while  the  associated  parameters 
are  summarized  in  Table  4. 


Table  4 


Parameters 


Clock  | 

1 

1 

1 

1 

b 

r  -13, 

[10  ] 

I  1 

1  C  1 

-1  S 

|  [10  /yr]  | 

1  1 

uw  d 

r  -13  i 
[10  ] 

1 

e 

1 

I  [d] 

1 

Csl  | 

1 

t— j 

o 

1  1 

1  -2.5  | 

-0.6 

1 

|  30 

Cs3  | 

-2.9 

1  0.9  | 

-0.5 

|  -30 

Cs5  | 

o 

00 

1  -0.2  | 

1.2 

1  o 

Cs6  | 

2.5 

1  -0.4  | 

1.1 

1  o 

Cs9  | 

1 

15.2 

1  0.6  | 

-2.4 

|  40 

1 
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The  coefficients  b,  c,  d,  e  are  those  found  in  equation  (4)  and  represent  respec¬ 
tively  : 

-  b:  the  mean  frequency  departure  at  x  =  0, 

-  c:  the  frequency  drift, 

-tod:  the  amplitude  of  the  sinusoidal  term, 

-  e:  the  initial  phase  (divided  by  uu)  at  x  =  0. 

In  the  lower  side  of  Fig.  6  and  Fig.  7,  the  temperature  and  the  absolute  humidity 
variations  for  both  IEN  and  ISPT  laboratories  have  also  been  reported. 

The  analysis  through  the  Allan  and  the  Hadamard  variances  was  performed  only  on 
the  Cs6  and  Cs9  IEN  standards  for  the  last  two  years.  The  statistical  data 
reported  here  below  are  extracted  from  reference  [ 1 3 ] . 

The  Hadamard  variance  was  used  to  remove  possible  linear  frequency  drifts,  but 
the  results  obtained  did  not  give  any  more  information  on  periodic  oscillations 
and  on  medium  term  instabilities.  As  a  consequence  from  the  data  of  both  the 
Allan  and  Hadamard  variances  it  was  suggested  that  frequency  drifts  were  negligi¬ 
ble  in  Cs6  and  Cs9  in  the  sample  time  considered  and  this  conclusion  was 
confirmed  through  the  Gauss-Newton  method  too. 

Cs9  Allan  variance  diagrams  referring  to  two  different  starting  times  for 
computations  are  shown  in  Fig.  8  where  seasonal  variations  appear  in  the  range 
25-250  days.  The  time  difference  between  the  start  points  in  the  two  graphs  is  90 
days  which  corresponds  to  AS0  close  to  -ff  /2  and  they  were  chosen  so  that  9a 
values  are  close  to  0  and  tt/2.  Note  that  in  this  figure  the  horizontal  axis 
reports  days,  one  day  being  approximately  1  degree.  The  minimum  and  maximum 
values  occur  for  sampling  times  equal  to  170  days.  Some  irregularities  may  be  as¬ 
cribed  to  the  small  number  of  data  available  for  the  averages.  In  effect  the  90 
dependence  of  cTy  ,as  obtained  from  many  graphs,  is  shown  in  Fig.  9  and  Fig.  10 
for  Cs9  and  Cs6  respectively.  It  is  strong  and  nearly  equal  for  X.  =  170  and  180 
days  whereas  it  decreases  for  X  =  160  days. 

A  rough  estimate  of  seasonal  variations  is  obtained  by  using  the  very  simple  ex¬ 
pressions  (8a)  and  (8b)  which  strictly  refer  to  a  deterministic  periodic  process. 
The  factor  =  (tV9*/  **4)* l<  *  in  the  non-stationary  case  =  TT  can  be 

easily  evaluated.  In  fact  /J<l turned  out  to  be  2.5x10  and  1.1x10  for 
Cs9  and  Cs6  respectively.  These  results  are  in  agreement  with  those  reported  in 
Table  4,  noting  that  uj  d  of  expression  (4). 

Some  conclusions  can  be  argued  about  the  cloks  analyzed.  Looking  at  Fig.  7,  where 
the  mean  frequency  departures  of  Csl  and  Cs3  of  ISPT  together  with  temperature 
and  humidity  diagrams  are  reported,  and  analyzing  the  coefficients  of  Table  4,  it 
is  suggested  that  the  amounts  of  their  frequency  variations  (cod  of  the  order  of 
1x10  13  peak  to  peak)  cannot  be  easily  related  to  the  humidity  variations  only, 
these  last  being  rather  small  whereas  the  temperature  varying  about  0.5°C. 
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On  the  other  hand  the  ISPT  clocks  synchronization  link  used  to  relate  them  to  UTC 
is  considerable  noisier  than  that  used  for  IEN  clocks  and  therefore  the  interpo¬ 
lation  results  have  a  low  confidence  level.  As  far  as  the  IEN  cesium  clocks  are 
concerned  (Fig.  6) ,  the  amounts  of  seasonal  fluctuations  range  from  2  to  5  times 
10  .  If  temperature  variation  effects  are  supposed  to  be  negligible  compared  to 
those  of  humidity,  we  can  ascribe  to  these  latter  the  observed  frequency  varia¬ 
tions.  Under  this  hypothesis  the  coefficients  of  fractional  frequency  variations 

—15  —3 . 

per  absolute  humidity  variations  (10  /g  .m  )  are  evaluated  to  be  24  for  Cs5,  22 
for  Cs6  and  -48  for  Cs  9.  These  results  are  in  agreement  with  the  evaluations 
reported  in  [2].  Moreover  a  delay  between  frequency  and  humidity  variations, 
ranging  from  some  days  to  one  month  and  depending  on  the  type  of  clock,  has  been 
pointed  out.  According  to  this  observation  it  seems  possible  to  find  out  a  suit¬ 
able  time,  for  each  clock,  in  order  to  compensate  the  seasonal  effect  and  to  im¬ 
prove  the  uniformity  of  the  IEN  time  scale. 

6.  Conclusions 

In  this  work,  the  presence  of  seasonal  effects  on  five  commercial  cesium  clocks 
has  been  detected  and  the  relative  amount  has  been  evaluated. 

Two  different  approaches  of  evaluation  have  been  experimented  obtaining  consis¬ 
tent  results  and  a  good  agreement  has  been  also  found  with  previous  investiga¬ 
tions.  At  this  point,  it  appears  possible  to  elaborate  suitable  criteria  of  a 
time  scale  seasonal  effect  compensation  and  to  maintain  UTC (IEN)  within  2  micro¬ 
seconds  . 
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CLOCKS  ROOM  TIME  AND  FREQUENCY  LABORATORY 


Fig.  1  -  Equipment  used  to  realize  UTC(IEN)  and  to  perform  clocks  intercomparisons 


Fig.  2  -  UTC(IEN)  vs.  UTC  from  BIPM/BIH  Annual  Report  and  circulars  D. 


Fig.  3  UTC(IEN)  vs.  UTC  compared  with  UTC(Cs6)  vs.  UTC.  The  mean  rate  of  Cs6 
is  subtracted. 
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Fig. 


4  -  IEN  clock  room  temperature  (T)  and  absolute  humidity  (Hu)  diagrams. 


Fig.  5  -  Normalized  Allan  and  Hadamard 

variances  for  periodic  modulation 
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Fig.  7  -  Normalized  frequency  departures  of  ISPT  cesium  clocks  vs.  UTC,  tempera 
ture  (T)  and  absolute  humidity  (Hu)  of  the  clocks  room. 


